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Abstract 

In this paper the authors describe a theoretical simple statistical modelling of relaxation process 
in metal-oxide semiconductor devices that governs its degradation. Basically, starting from an 
initial state where a given number of traps are occupied, the dynamics of the relaxation process is 
measured calculating the density of occupied traps and its fluctuations (second moment) as function 
of time. Our theoretical results show a universal logarithmic law for the density of occupied traps 



(n(t)) ~ <p(T,Ep) (A + B hit), i.e., the degradation is logarithmic and its amplitude depends 
on the temperature and Fermi Level of device. Our approach reduces the work to the averages 
determined by simple binomial sums that are corroborated by our Monte Carlo simulations and 
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by experimental results from literature 



which bear in mind enlightening elucidations about the 



physics of degradation of semiconductor devices of our modern life. 
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The understanding of the physics of semiconductor devices [l| has never been so impor- 
tant, since silicon-based integrated circuits are facing increasing reliability and scaling issues. 



Quantum Computing [16|, DNA Computing [15|, and many other alternatives are arising as 



possible substitutes to the technology of the silicon-based computation, but these seem to 
be very distant from the reality of our day-by-day. 

Hence, the understanding of the reliability effects in semiconductors is of paramount 
importance for the Microelectronics Industry, that is more and more dependent of new 
results from an interdisciplinary Physics and their ramifications. 

Complementary metal-oxide semiconductors (or simply CMOS devices) have an impor- 
tant role in the development of the information and electronic industry. A scheme of this 
device can be represented basically by figure [2J As can be observed, it is composed by a thin 
metal plate, followed by an insulating layer (e.g., SiC^), and finally a semiconductor layer 
(Silicon-Si). The working of device supposes two voltages: 

1. The gate-source voltage (Vgs), which controls the Fermi- level of charge carriers in the 
semiconductor, i.e., the electrons are attracted with larger or smaller intensity to the 
interface between the oxide and semiconductor according to the magnitude of Vqs] 

2. The drain source voltage (Vds), responsible to move the electrons to compose the 
current (drain current). The uniformity of this current can be affected by charge 
traps found close to the semiconductor-insulator interface. In this context the rules 
of capture and emission of charge carriers by traps in semiconductor devices lead to 
irregular signals (noise) in the current which can be observed in figure [TJ 

In this figure we observe how the current of device is affected by a single trap. The fluctu- 
ations experimentally observed are caused by the effect of the superposition of a number of 
traps, under different conditions and moreover, we can consider fluctuations from a sample 
of devices for completely describing this phenomena (see for example |4|]). 



Such fluctuations are called random telegraph signals 131. In other contexts, as for example 

n 

Ni/NiO/Co junctions [10|], such fluctuations are even able to govern magnetoresistance and 
;heir study also will be very important in the context of nanostructures (see for example 



(12|). 



The capture and emission of charge carriers by the traps is described as a simple Poisson 
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FIG. 1: A simple representation of a random telegraph signal caused by sucessive captures and 
emissions of charge carriers by a trap in the current of a CMOS transistor. 

process governed by rates r c and r e , where the capture occurs with probability p(0 — > l)dt = 
r~ x dt and and emission p{l — > 0)dt = r~ x dt 

Such rates r c and r e can be defined as the time average in state 1 and state respectively: 
r c = (t)i — Jq° t Pi(t)dt and r e = (t) = J °° t Po(t)dt, where Pi(t) is the probability 
of permanency in state 1 and Po(t) the respective amount for state 0. Naturally Pi{t) = 
l/r c ■ e" t/rc and P (t) = l/r e ■ e^ 1 ^ . 

If the number of trapped charge carriers increases over time, a decrease of the current 
may be observed. This is an aging effect usually called bias temperature instability (n- or 
p-bti), since it depends on bias (Fermi level) and temperature, as discussed below. 

In this context the degradation of a MOS transistor can be measured as the number of 
occupied traps and the dynamics of this occupation must be better understood. In this paper 
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FIG. 2: A scheme of a complementary metal oxide semiconductor device (CMOS transistor), 
composed by a fine metal plate, an oxide (Si02), and a semiconductor (Si). A voltage Vgs (gate 
source) has a role of attracting the electrons to the top of semiconductor near of the oxide. Other 
voltage Vds moves the charge carriers generating a drain current. The charge carriers under a 
Fermi Level can be captured or emitted come back to the current along time evolving according 
to Fermi-Dirac statistics. These successive captures and emissions generate the known random 
telegraph signals. 

we aim at a theoretical analysis to describe the density of occupied traps in a semiconductor 
device and so to understand how a characteristic degradation process occurs in these devices 
and other similar devices. 

So, first of all, we need to calculate the probability of a particular trap with constants 
r c and r e starting from state (empty) and after a elapsed time t it is in this same state, 
which we denote poo(t). This probability can be calculated observing that ([lj], jsj): 

P 01 (t + dt) = P„i(t)p(l -> 1) + ^oo(*)p(0 -> 1) 
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where p(0 — > 1) = dt/r c and p(l — > 1) = 1 — p(l — >■ 0) = 1 — dt/r e and also Poo(t) = 
1 — Poi(t). This leads to a simple differencial equation: dPo\{t)/dt = r" 1 — (r" 1 + r ( 7 1 )P i(^)- 
If P o(0) = 1, its solution is P i(t) = r e (r e +r c ) _1 [1 - exp(-t/r e<? )], where l/r e(? = l/r e +l/r c . 
Similar results can be performed leading to -Pn(t) = — \j e + r c exp(— t/r eq )}. Unless 

T e + T c 

fluctuations on the current amplitude, whose average A depends on other microscopic factors 
of the device, this probability corresponds to the autocorrelation of the system A(t) = 
(a(0)a(t)), where a{t) corresponds to the state of the trap (occupied a = or empty 
o" = l). In frequency domain this exponential decay of autocorrelation for one trap is 
described by Lorentzians, since the power spectrum density, i.e., the fourier transform of the 
autocorrelation is 



/oo 
e -^f«A(t)dt 
-oo 



4A 2 



^ + r c )(4 + 2vr/2) 

The known 1/ / noise results from a sum of these Lorentzians (a contribution o 
Taps in device). For more details about the origin oil/ f noise, see for example 
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So, coming back to relaxation phenomena, and starting from n(0) = (all traps empty), 
we can calculate the average density of occupied traps at time t 

(n(t)) = {Et r o^(t) / 

(1) 

where Pr(fc|n(0), t) is the probability of just k traps are occupied at time t, with k = O...N tr - 
But the traps have different constants r e and r c and from that, we write 

k N tr 

Pr(*|n(0),f) = X;n p oi(rj*\rj*>;t) J] P m {r^\ r<*>; t) 

C k i=l i=k+l 

with Ck denoting every subset {di, di, df-} from {1, 2, , n}. But 
statistically independent and identically distributed, and we have : 

k TV _fc 

Pr(%(0),t) = Ec, P oi(r c ,r e ;t) • P m (r C} r e] t) " 



(l)Poi(T c ,T e ;t) ■ P o(r c ,r e ;t) 



■Ntr-k 
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where 7 = f f ■ f(T c )g(T e )dr c dT e , where f(r c ) and g{j e ) are probability densities of time 
constant times of capture and emission. Microscopically these quantities can be better 
understood. Actually, r c and r e have a thermal and a gate voltage dependence. Some more 



detailed approaches use quantum two-dimensional calculations of these quantities 



9|. Here 



we use a known simplification proposed by Kirton and Uren [1] , where r c and r e are random 
variables that follow the form: r c = 10 p (l + exp(— q)) and r e = 1(F(1 + exp(g)), where 
P e [j5 min ,p max ] and q = (E t - E F )/k B T 6 [E v - E F )/k B T, (E c - E F )/k B T] are randomly 
distributed according to respectively a uniform and a u-shape distribution. 

There is no much information about these density of states of the traps in literature, but 
Wong and Cheng (| show that for 3 different prepared gate oxides it follows a u-shape form. 



Naturally, we must observe that r, 



<<i 



r c T e /(r c + r e ) = 10 p corresponds to an uniform 
distribution of time constants (r eq ) in a log scale, as expected. Here, E t and Ep are respec- 
tively the energy of observed trap and Fermi Level of system that is directly proportional 
to Vqs applied in device. For our purpose, it is more interesting to switch our average: 




') f(T c )g(r e )dT c dT e ->■ 

max Pmin ) 

where fl(E t ) is the density of states of the traps in the interface. 

Coming back to equation [TJ after some straithforward calculations we have 



n(t)} = Ek=oHl) P oi(P,Eut) ■ P Q i(p,E t ;t) 
= N tr P 01 (p,E t ;t) 

- ISI trJ Ev 1+e -(E t -E F )/k BT a^f 

-, l - T f Pmax dp [1 - exp(-10- p t)l 

(Pmax-Pmin) J Pmin r 1 ^ K Ji 



Ntr 



The second integral must be better worked out. Making a suitable change variable p = 
— log 10 (j), dp = — In -1 10^ and we have the temperature dependence separated of time 
dependence via two integrals: 

r E ° dE t n(E t ) 

l E 1 _|_ e -(E t -E F )/k B T 



(n(t)> = N t 



tr 



In" 1 10 



Pn 



1Q-P max t 



du- 



1) 



10-Pmint 



U 



(2) 
(3) 



which can be analyzed numerically. A particular case is when f(E t ) is uniform, and in this 



case we 



have j 1 



E c dE t g(gt) 

E v l +e -( E t- E F)/ k B T 



k B T 



(E c —E v ) 



In 



e (Ec-Ev)/k B T +e (E F -Ey)/k B T 
l+e^F ~Ec)/k B T 



. HE F ={E v +E c )/2, 
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i.e., it is exactly in the middle of band gap, this integral is numerically equal 1/2 and there 
is no temperature dependence, i.e., 



Ht)) umform = 7u~ / du K - } (4) 



I uniform 9/ _ \ / "" 

A simple particular case is if we observe the evolution of the occupation probability of a 
single trap, with time constants r c and r e , is numerically equal to lim^-^oo ( n (t)) uniform/ '^tr 
when all traps have the same r e and r c , which is given by: 



1 = — ■ 



Pr(^(t) = 1) = Ytp- (5) 

where /3 = r c /r e is a important ratio considered in this context. Naturally Pr(<7j(t) = 1) — > 
j^p = - T ^_ r when t — > oo. This leads to a simple but important conclusion: if /3 > 1 
(time capture is greater than time emission) we have Pr(cjj = 1) < 1/2 and otherwise - i.e., 
f3 < l(time emission greater than time capture), Pr(crj = 1) > 1/2 when t — > oo. 



But what is the behavior of (n(t)) in a realistic case(i.e., when r c and r e are randomly 
distributed)? In this case we solve numerically the exponential integral from El We adopt 
usual values found in the literature for this problem (p m i n = and p max = 7) what means a 
frequency ranging from 1 to 10 7 Hz. 



The continuous curve in figure [3] shows the time evolving (n(i)) theoretically obtained 
(i.e. numerical integration of equation |2]). 

The points corresponds to our MC simulations. For these MC simulations, we start from a 
given number of empty traps. For each time, each empty trap % = 1, N tr becomes occupied 
with probability pi(0 — » 1) = 10~ Pt (l + e _9l ) _1 and similarly an occupied trap becomes 
empty again with probability Pi(l — >■ 0) = 10 _Pl (l + e 91 )^ 1 , where Pi is uniformly drawn in 
bmm,Pmax], and here the same values (p m \ n = and p max = 7) were used. Similarly is 
uniformly drawn in [E V ,E C ]. For real devices E t should change from ~ 0.2 eV to ~ 1.0 eV. 
First of all, Ep corresponds to middle of band-gap (Ep ~ 0.6 eV), ks = 8.617385 x 10~ 5 eV 
K _1 and T = 300 K, which leads to qi G [—15.473,15.473]. We can observe a excellent 
agreement between the MC simulations and our theoretical equations, showing that in semi- 
log plot in figure |3] the relaxation dynamics follows a logarithm law: 

Ht)) uniform ~ A + B\ogt 
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FIG. 3: Time evolving of (n(t)) obtained directly by numerical integration of equation [5J The points 
corresponds to MC simulations performed in same conditions. 

where we find A = 31.80(15) and B = 144.20(3). The uncertainties were obtained, using 
error bars obtained from 16 independent runs of the program. We test other temperature 
values but we did not find difference as expected when Ep is in the middle of band gap in 
this case where a uniform density of states is considered. 

Extending our results, we can compute the second moment: 



N tr 

; n (t) 2 ) = ^£; 2 Pr(%(0),£) 

k=0 



•2 



which yields (n(t) 2 } = N tr Pr(k\n(0),t) + N tr (N tr - l)Pr(fc|ra(0),t) « (n(t)) + (n(t)) 



So we have {n{tf) « 10 , f*" dJ-^ + 10jV ^ ' 

\ \ / / 2(p max -p min ) J 10 'mini u 4(p max -Pmin) 

which leads to 



10 -pmaxj ( e -"_lj 
JlO-Pmint dU U 
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FIG. 4: A polynomial fit of a u-shape (Reoxidized Nitrided-oxide) . The experimental curves were 
extracted from (|. The inside plot corresponds to the other u-shape densities that are very similar. 



{n(t) 2 ) ~ A(l + A) + B(l + 2A) logt + B 2 log' t 

with A and B exactly as reported before. 

Looking at the temperature dependence, we must analyze more realistic densities of 
states. After a detailed scanning of the plots for the densities of states, for the 3 prepared gate 
oxides (TCE Oxide, Reoxidized Nitrided-oxide, and Nitrided-oxide) found in the reference 
6], a fitting by a eighty-degree polynomial, here described by Vt(E t ) = Ylk=o^ k ^t were 
performed (see {2] for a more detailed discution of this part). This excelent fit can be seen 
in figure HI 

Using these u-shaped densities of states or even their polynomial fit, we can calculate the 
temperature dependence: 
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x 1 v^ 8 n f Ec E t d Et 
<pF, E F ) = - (3 k ^ 1 + J E J F)/kBT (6) 

with Z = Yl=^k{k + - E* +1 ) where then (^(t)} ~ ^(T,£ F )(A + Blogt), 

showing that temperature and Fermi level dependence are independent of time, i.e., the 
relaxation depends logarithmically on time for each fixed pair (T, Ep). 

Using Simpson numerical integration, we calculate <p{T) and now the time evolving of 



(n(t)} is calculated according to equation |2j By motivations from experimental results, 



initially we study (n(t)) for a Fermi level set on Ep = Ec = leV using the u-shape ob- 



tained from Reoxidized Nitrided-oxide. The figure [5] shows (n(t)} for different tempera- 
tures. The inside plot in this same figure shows the amplification factor of temperature 
A(t) = ip(T)/(p (323), since T = 323K (50°C) is the minimum temperature used in our 
calculations, for the 3 different oxides extracted from 6|. We show the linear universal 
behavior, described by relation 

A(T) ~ 0.18 T 

It is also interesting to study the dependence on Fermi level. The Fermi level may be 
varied by changing the gate bias of the device as experimentally explored in [5j . We studied 
the dependence of tp{T) as function of T for different values of Ep. Our results show that the 
curves showing the evolution of relaxation as function of time can be collapsed in a single 
one if multiplied by a suitable constant. This constant is (p(T,Ep) (see figure E]). Hence, 
our results are in agreement to the experimental findings for the Fermi-level dependence of 
the relaxation. 

We did also study the dependence on the initial density of occupied traps. In this case, if 
the density of initially occupied traps is lower than the equilibrium value, in the relaxation 
process, the number of occupied traps increases logarithmically. On the other hand, if 
the density of initially occupied traps is higher than the equilibrium value, the density 
of occupied states decreases logarithmically in the relaxation process. MC simulation did 
confirm this behavior. MC simulations were performed starting from different initial density 
of occupied traps p = 0, 0.1, 0.2, 0.3, 1.0, where the equilibrium value is p = 0.5. We 
can observe this behavior in figure [71 

In summary our results corroborate the experimentally observed logarithmic relaxation 
of the density of the occupied traps via MC simulations and from theoretical analysis in 
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FIG. 5: Density of occupied traps for different temperatures (323, 373, 423, 448, 473 K) using a 
real u-shape,we observe effects on logarithmic law exactly as observed via experimental results 
obtained in 5]. The inside plot shows the amplification <p(T)/<p(323) as function of temperature 
for the 3 different oxides showing a agreement among them indicating a linear universal behavior 
ip(T)/ip(323) ~ 0.18T. 

complementary metal-oxide semiconductors governed by Fermi-Dirac-Shockley-Read Statis- 
tics [8]. Our results also corroborate the experimentally observed temperature dependence, 
which shows that the relaxation as a function of time at different temperatures may be 
collapsed into a single curve using a suitable scaling factor. This behavior is experimentally 



observed e.g. in The scaling factor is the function ip{T,Ep) of equation El 
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